Digital Quantum Simulation of the Statistical Mechanics of a Frustrated Magnet 
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Many interesting problems in physics, chemistry, and computer science are equivalent to problems 
of interacting spins [T]. However, most of these problems require computational resources that are 
out of reach by classical computers. A promising solution to overcome this challenge is to exploit 
the laws of quantum mechanics to perform simulation [2]. Several "analog" quantum simulations 
of interacting spin systems have been realized experimentally [3-8 . However, relying on adiabatic 
techniques, these simulations are limited to preparing ground states only. Here we report the first 
experimental results on a "digital" quantum simulation on thermal states; we simulated a three-spin 
frustrated magnet, a building block of spin ice, with an NMR quantum information processor, and 
we are able to explore the phase diagram of the system at any simulated temperature and external 
field. These results serve as a guide for identifying the challenges for performing quantum simulation 
on physical systems at finite temperatures, and pave the way towards large scale experimental 
simulations of open quantum systems in condensed matter physics and chemistry. 



The most challenging aspect of many-body simulation 
is that the memory and temporal resources often scale 
exponentially, rendering many problems of interest in- 
tractable by all known classical methods [2 . A promis- 
ing solution is quantum simulation^ in which a quantum 
system acts as a processor to simulate another physical 
system (quantum or classical). There are two classes of 
quantum simulation: "analog" simulators are typically 
engineered to simulate a particular class of Hamiltoni- 
ans [9] and to find ground states of non-trivial Hamil- 
tonians adiabatically, whereas "digital" simulators rely 
on universal quantum information processors (QIPs), ca- 
pable of implementing a universal set of quantum gate 
operations [10] to simulate not only ground states, but 
also thermal states, and even time-evolving states. 

Simulations of interacting spin systems are of partic- 
ular importance to many applications, such as modeling 
magnetism |TT], solving optimization problems [12], and 
restoring digital image [13]. Furthermore, understand- 
ing the properties of the spin models also offers insights 
to the computational complexity theory [M]. For exam- 
ple, the ground-state problem of the Ising spin model is 
known to be an N P-complete problem; this implies that if 
an efficient algorithm for solving the ground-state prob- 
lem of the Ising model exists, then it can solve all other 
problems in the class of NP. This matter is related to the 
question whether P equals NP, and is a major unsolved 
problem in computer science. 

In a series of recent experiments [3-6 , based on the 
adiabatic methods, progresses of quantum simulation of 
various spin systems have been achieved in various phys- 
ical systems. These experiments, however, suffers from 
two limitations: (a) they are limited to studying the 
ground-state properties only, and (b) the energy gaps 



along the adiabatic paths must be large enough to avoid 
excitations from the ground states. In general, the en- 
ergy gaps cannot be pre-determined efficiently, and are 
non-controllable. Therefore, the advantage of the adia- 
batic methods over classical methods is not guaranteed 
for all cases. 

On the other hand, at finite temperatures, all of the 
thermodynamics of spin systems can be obtained by de- 
termining the partition function which, instead of NP, 
falls into a different complexity class called sharp- P, or 
#P. However, if an efficient algorithm for evaluating par- 
tition functions exists, then the ground-state properties 
of the corresponding spin systems can also be determined 
efficiently. Therefore, the problem of determining parti- 
tion functions is at least as hard as the NP-problems, or 
simply called NP-hard. 

Practically, partition functions cannot be computed ef- 
ficiently, except for some simple cases such as 1-D spin 
chains. For classical spins, the classical Metropolis al- 
gorithm provide a means for generating the Gibbs dis- 
tributions, through the construction of Markov chains 
with Monte-Carlo methods. For quantum systems, the 
quantum generalization of the Metropolis algorithm has 
been achieved JSl [16]. However, Markov-chain based 
methods, similar to the adiabatic methods, are limited to 
the cases where the Markov-matrix gaps cannot be too 
small to achieve convergence. Particularly, for frustrated 
spin systems. Metropolis sampling can result in ensem- 
bles trapped in local minima. In these cases, methods for 
direct encoding the Gibbs distribution into the states of 
the qubits would be more efficient. This is the key issue 
that motivates this experimental work. 

By using a digital simulator, our system can explore 
the full phase diagram of the thermal state of the sys- 
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tern, as a function of temperature and magnetic field. 
All ranges are in principle accessible due to the digital 
nature of the simulation. Unlike other quantum phase 
transitions explored in other quantum simulation experi- 
ments ^-i6j , the underlying Hamiltonian of the simulator 
does not need to be restricted to certain types of inter- 
actions to be able to simulate the target system. 

In this letter, we report the first digital quantum sim- 
ulation of the finite-temperature properties a classical 
three-spin frustrated magnet, a building block of spin ice, 
(see Fig. [T]), using a four-qubit quantum register based 
on NMR. The reason for simulating the frustrated mag- 
net is that it exhibits a rich phase diagram of the total 
magnetization as a function of temperature and magnetic 
field. This allows us to experimentally probe various dis- 
tinct features of this system. On the other hand, the 
phenomenon of geometric frustration is an interesting 
topic in condensed matter physics. For example, mate- 
rials, such as water ice, exhibiting geometric frustration 
cannot be completely frozen; the motion at the molecu- 
lar scale continues even at absolute zero. Recently, the 
same three-spin frustrated magnet at zero temperature 
has been simulated by trapped ions [5[ E]. We aim to 
make progress along this direction by extending the quan- 
tum simulation of the frustrated magnet to finite temper- 
atures. 

In our simulation, instead of a mixed state, the im- 
plemented algorithm prepares a coherent encoding of 
a classical thermal state (GETS) on a quantum regis- 
ter [HI HI], 
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which is a pure state (a pseudopure state in the NMR 
experiment) with amplitudes ^/e~^^k/Z equal to the 
square roots of the corresponding thermal state Gibbs 
distribution associated with the eigenstate of the 
Hamiltonian H. Here P = 1/T {ks = 1), and Z = 
Tr(e-^^) is the partition function. The GETS, there- 
fore, contains all of the information about the thermal 
density matrix 
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of the system. In fact, the thermal density matrix pth 
can be directly obtained from the GETS state \^ ^) by 
artificially "decohering" the off-diagonal elements of the 
density matrix \^ ^)(^ ^\ constructed from the GETS. 

In this method [18 , the number of quantum gates 
needed to prepare the GETS is linear in the number of 
spins for ID cases and sub-exponential for 2D cases, but 
is still exponential in general for NP-problems. Nonethe- 
less, the efficiency of this algorithm is independent of the 
simulated temperature, and not limited by the small-gap 
problem encountered in the Markov-Ghain Monte Garlo 
algorithms. This makes it advantageous for simulating 
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FIG. 1: Theoretical descriptions of the frustrated magnet, 
(a) All possible configurations of a three-spin frustrated mag- 
net at zero temperature and zero magnetic field. There is a 
six- fold degeneracy in the ground state, leading to a non-zero 
entropy, (b) Theoretical phase diagram. The units of the axes 
are ku and J for temperature and external field. The dashed 
line parallel h- axis denotes T — 1/11, which corresponds to 
the case of the experimental data taken in fig. 3. (c) Quan- 
tum circuit diagram for preparing and measuring the GETS 
1^^) defined in Eq. [ijfrom the initial state |0000). The gate 
sequence for preparing the GETS on the lower three qubits 
is determined by the method in Ref. ^8J. An explicit way 
for constructing the [/fc -gates is given in the supplementary 
materials. The top qubit ^0 serves as a probe for measuring 
the physical observables. 



the low-temperature properties of frustrated spin sys- 
tems. Furthermore, although this algorithms can at most 
yield a quadratic speedup for simulating the most general 
thermal states [18 , the subclass of the GETS which can 
be created efficiently on a quantum computer can serve 
as a "heat bath" [19 for the simulation of the dynamics 
of open quantum systems, which could give an exponen- 
tial advantage [20 . Our goal is to investigate how well 
such a GETS can be prepared in the laboratory subject 
to the existing experimental constraints. 

In this experiment, three qubits encode the GETS of a 
triangle plaquette of Ising spins with equal couplings J, 
with temperature T and global magnetic field h as vari- 
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ables. A fourth ancilla qubit is used to probe the physical 
properties of the GETS by measuring the set of diagonal 
Pauli operators, so that quantities such as the total mag- 
netizations and spin correlations can be extracted. These 
measurements are sufficient for determining the partition 
function from which any thermodynamic quantity of 
interest, e.g. entropy can be calculated. 

To be more specific, the Hamiltonian of the frustrated 
magnet is defined by 

U = J (Z1Z2 + Z2Z3 + Z1Z3) + (Zi + Z2 + Z3) , (3) 

where Zi = is the z Pauli matrix of the spin i. For 
J>0, the coupling is antiferromagnetic, where the spins 
tend to minimize the energy by pointing in opposite di- 
rections. The external field /i, however, tends to force 
the spins to align. The effect of finite temperature T 
is to wash out both tendencies. The competitions be- 
tween these factors give rise to a phase diagram with a 
rich structure, as shown in Fig. ^p. For example, near 
some critical values of the external fields h = —2 J, and 
2 J, there are crossover points where the configuration 
of the spins, and hence the total magnetization, change 
abruptly. Near T = and = 0, the ground state is 
fully frustrated with a six-fold degeneracy, as illustrated 
m Fig. [1^. This means that, unlike ordinary materials, 
the entropy (and hence heat capacity) of the frustrated 
spin system is non-zero at T = 0. 

Our strategy for the study is as follows: for any 
given value of the temperature T and magnetic field /i, 
the GETS 1^/3) can be prepared with a quantum cir- 
cuit of constant depth, as shown in Fig. [l]^. Here we 
choose J as the unit for h. The three lower qubit s, 
qi,q2 and gs, initialized into the |000) state, are cho- 
sen as register qubits to encode the GETS. We choose 
the phase kick-back method [21 to extract information 
about the GETS by introducing a fourth qubit qo as a 
probe qubit, shown as the top qubit in Fig. [l]^. The 
probe qubit is then prepared in a superposition state 
(|0) + |l))/>/2 and a controlled- /7m gate operation is ap- 
plied to the joint probe-GETS system to measure observ- 
ables (Um) = (^/3|^m|^/3) on the top "probe" qubit. 
Here (Um) is proportional to the coherent element in the 
reduced density matrix of the probe qubit. 
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through tracing over the register qubits and can be ex- 
tracted from the NMR signal of the probe qubit. By 
measuring the set of operators 

Um = {^1j^2j ^3, ^1^2, ^2^3, ^l^Sj Z1Z2Z2,} , (5) 

the full thermal state density matrix pth can be recon- 
structed from the GETS. 

For the NMR implementation, we choose as a sam- 
ple of the ^^G- labelled trans-crotonic acid dissolved in 



d6-acetone, which forms a seven-qubit register; the four 
qubits in this experiment corresponds to the four carbon 
spins, and the other three nuclear spins are not directly 
involved after the preparation of the pseudopure state. 
The experiments were carried out on a Bruker DRX 700 
MHz spectrometer. The structure of the molecule and 
the Hamiltonian parameters of the seven spin qubits are 
shown in Fig. [2^, where the NMR Hamiltonian of this 
system is given by 



Hnmr = -tt ^ i^icrt + TT ^ Jkicrl.af/2 
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where Ui denotes the chemical shift of spin i, and Jki 
denotes the coupling strength between spins k and 

In the experiment, we exploit standard Isech and 
Hermite-shaped pulses to implement single-spin opera- 
tions for the nuclei M and G1-G4, and numerically op- 
timized GRAPE pulses |24] for manipulating Hi 
and H2 (for initial pseudopure state preparation only). 
A custom-built software compiler generates pulse se- 
quences, including refocussing pulses, that are optimized 
for the highest unitary fidelity [24 . Furthermore, the 
radio- frequency (r.f.) spin selection techniques are ex- 
ploited to improve the linewidth, and therefore the co- 
herence, of the ensemble qubits [22l [25]. The effect of 
pulse errors due to r.f. inhomogeneity is reduced by a 
spatial selection of molecules in a region of high r.f. ho- 
mogeneity. A labelled pseudo-pure state is prepared of 
the form ps = OOa^OOOO using methods described in Ref. 
|22J, where = |0)(0| and the order of qubits is as fol- 
lows: M, Hi, H2, Gi, G2, G3, G4. Note that we employ 
the deviation density matrix formalism [27 . 

The four carbon spins, initialized in the state 0000, 
are used to prepare and measure the GETS, where Gi 
is the probe qubit, and G2 - G4 are the register qubits 
for simulating the frustrated magnet. The GETS 1^/3) is 
prepared by the pulse sequence shown in Fig. ^jp. The 
NMR signal of Gi is acquired after the controlled- /7m 
gate is applied. The controlled- /7m is implemented by 
combining phase-fiip and SWAP gate operations, and can 
be further be decomposed into nearest-neighbor coupling 
evolutions and single spin rotations. Fig. ^ illustrates 
the sequence for the observable Z1Z2Z3. 

As indicated in Eq. Q, (/7m) is encoded in the co- 
herent part of the probe qubit (Gi) state. In the spectra 
of probe qubit (see supplementary material), the coher- 
ence is distributed among 2^ = 64 peaks, each of which 
corresponds to a particular eigenstate of the remaining 6 
qubits M, Hi, H2, G2-G4. The intensities of these peaks 
are obtained by a precise spectral fitting procedure [28] . 
In preparing and measuring the GETS, no computational 
operations are performed on the proton spins M, Hi, and 
H2. The numerical simulations of the experiment take 
into account the effect of Ti relaxation process of the 
proton spins, which is experimentally measured though 
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FIG. 2: Experimental protocol, (a) Hamiltonian parameters 
for the nuclear spins in carbon-13 labelled trans-crotonic acid, 
with the structure shown as the inset. The chemical shifts and 
J-coupling constants (in Hz) are listed on and above the diag- 
onal in the table, respectively. The longitudinal and transver- 
sal relaxation times Ti and T2 measured by standard inver- 
sion recovery and Hahn echo pulse sequences are listed at 
right. The chemical shifts are given with respect to reference 
frequencies of 700.13 MHz (protons) and 176.05 MHz (car- 
bons). The molecule provides seven qubits since the methyl 
group can be treated as a single qubit using a gradient-based 
subspace selection [22 . (b,c) Pulse sequences for preparing 
the GETS and measuring Z1Z2Z3 via the probe qubit, re- 
spectively, where refocusing pulses are not shown. The three 
carbons G2, G3 and G4 act as the GETS register qubits 1, 
2 and 3, and Gi acts as the probe qubit. The rotation an- 
gles are determined by the angles in Fig. [l]are listed in the 
supplementary material. 
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FIG. 3: Experimental measured magnetization and correla- 
tions, (a) Magnetization Zi + Z2 + Z3 as a function of mag- 
netic field h at low temperature, T — 1/11. The experimen- 
tal data ("o") is plotted together with numerical simulation 
results ("x") that include effects of carbon T2 and proton 
Ti. The theoretical result is shown as the solid curve. The 
points labeled "Experiment (R)" are obtained from the ex- 
perimental data by using a simple decoherence model to par- 
tially remove the effects of decoherence with no free param- 
eters (described in the text). The sharp change show the 
phase transitions, and the regions around the critical points 
h = —2, 0, 2 are enlarged as figures (b-d), respectively. (e) 
Surface plots for the total magnetization Zi + Z2 + Z3, corre- 
lations Z1Z2 -\- Z2Z3 + Z1Z3, and Z1Z2Z3 in theory (top row) 
and measured in experiment (bottom row). 



the decay of their initial state OOcr^llMll, where 11 de- 
notes the identity operator (see supplementary material 
for details). 

Since = — 1 with 1 = |1)(1|, the 64 peaks are 
divided into two antiphase multiplets corresponding to 
the two eigenstates of H2 (see supplementary materials) . 
We may, for example, choose the group marked by the 
H2 state 0. By adding the intensities of the 8 peaks 
marked by the state 0m^Hi^H2^ (Um) is obtained, taking 
into account proper normalization relative to the initial 
pseudopure state ps. 

The experimental results are summarized as follows: 
the diagonal elements of the density matrix constructed 



by the GETS are determined by measuring the full set 
of diagonal Pauli operators [see Eq. ([5|] for a range of 
simulated temperatures T and external fields h. The ex- 
perimental results are shown in Fig. |3] In FigjS^, for 
a low temperature (/3 = 11), the total magnetization 
Zi + Z2 H- Zs of the frustrated magnet is probed for a 
range of the simulated field h. The raw experimental 
results are in good agreement with the numerical simula- 
tion which takes into the account the decoherence effects. 
By rescaling the total magnetization by a constant fac- 
tor, which is equivalent to removing the isotropic errors 
^29j (see also supplementary material), the rescaled re- 
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suits agree much better with the theoretical predictions. 
In any case, it is clear that the magnetization changes in 
steps when the simulated magnetic field is varied from a 
large negative value {h = —5) to a large positive value 
{h = 5). The critical points for the crossovers are lo- 
cated ai h = ±2, and 0, in agreement with the theoret- 
ically prediction shown in Fig. [l^. This is the result of 
the competition between the antiferromagnetic couplings 
and the external field applied to the frustrated magnet. 

Apart from the total magnetization, the other correla- 
tion functions are probed systematically for a range of the 
simulated temperatures and external fields. The results 
are shown in Fig. [3^. From these data, we can construct 
the thermal state density matrix of the frustrated magnet 
(subject to the normalization condition Ti (pth) = !)• 

Pth = ^t^^aiZi ^^bjkZjZk + CZ1Z2Z3 , (7) 

i j<k 

where ai= (Zi) /8, bjk= (ZjZk) /8, and c= (Z1Z2Z3) /8. 
We ignore the imaginary parts of the elements, which 
are zero in theory, and less than 11% in the experimen- 
tal data. With complete knowledge of the thermal den- 
sity matrix pth^ we can determine all of the macroscopic 
thermodynamic observables for an ensemble of frustrated 
magnets. In this study, we are particularly interested in 
investigating a non-linear quantity, namely the entropy, 

S = -TT{pth\npth} , (8) 

and gauge how sensitive it is to experimental errors in 
{Um). 

Fig. [4^ shows the experimental results for the en- 
tropy 5 as a function of the simulated magnetic field h in 
the low temperature regime (/3 = 11). One finds that the 
sharp changes of S around h = ±2 and correspond to 
the crossovers, which have been observed in measurement 
of the total magnetization. Compared with the region of 
\h\ > 2, the large values of S in the region \h\ < 2, 
especially around = 0, indicates the preference of anti- 
ferromagnetism which causes the frustration for the mag- 
net. For the outer region \h\ > 2, where the external field 
should be strong enough to polarize the magnet, the the- 
oretical predictions of S should be zero. However, the 
experiment results, including the rescaled results, show 
non-zero values. This is due to the anisotropy of the mea- 
surement results taken from the three nuclear spins. The 
surface plots in Figj4]3-d show the entropy as a function 
of h and /3 from theory (b), experiment (c), simulation 
(d), and the rescaled experimental results that partially 
remove decoherence effects (e). We see that the role of 
the temperature is to "wash out" the competition be- 
tween the antiferromagnetic coupling and the external 
field. This is indicated by a transition near some value 
of the temperature (/3 = 3), beyond which the variation 
of the simulated external field h no longer causes sharp 
crossovers. 
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FIG. 4: (a) Entropy 5* as a function of magnetic field h at low 
temperature, /3 = 11. The experimental data ("o") is plotted 
together with numerical simulation results ("x") that include 
effects of carbon T2 and proton Ti. The theoretical result is 
shown as the solid curve. The sharp changes of S around 
h = ±2 and indicate the phase transitions. The points la- 
beled "Experiment (R)" are obtained from the experimental 
data by using a simple decoherence model to partially re- 
move the effects of decoherence with no free parameters (see 
SI), (b-e) Surface plots of entropy as a function of h and /3 
from theory (b), experiment (c) and simulation (d). Modified 
experimental results that partially remove decoherence effects 
are shown in (e). 



To summarize, the imperfection of the experimental re- 
sults is mainly due to decoherence effects; the duration of 
the experiment ranges from 0.35-0.76 s, (see supplemen- 
tary material for details), which is comparable with the 
T2 times of the nuclear spins. The simulated data, which 
take into account the nuclear T2 decay, closely match 
the experimental results. To a much lesser extent, inho- 
mogeneities of the applied magnetic fields and imperfect 
pulses also contribute to the error. Partial recovery of 
the density matrix is possible for isotropic errors (see 
supplementary material), this allows us to rescale the 
magnetization to improve the agreement with the theo- 
retical values (see Fig. [3|. The anisotropic part of the 
error, however, cannot be removed, and significantly af- 
fects the experimental entropy results (see Fig. [4|. 

The phase kick-back method was chosen for readout to 
exploit the fact that all J-couplings are well resolved for 
spin Ci. For spins C2 and C3, there is non-negligible 
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overlap between certain spectral peaks, which makes 
direct state tomography unreliable for certain states. 
Nonetheless, for certain simple states like |1111) at at 
/3 = 11 and = 5, state tomography through direct 
readout of all three spins can be performed reliably and 
compared to the phase kick-back results. The state fi- 
delity measured in this fashion was close to the phase 
kick-back result, only differing by 0.4% (see Fig. [g] in 
supplementary information). 

In conclusion, we employed a quantum information 
processor using nuclear spins to perform a digital quan- 
tum simulation of a geometrically frustrated magnet sub- 
ject to a simulated magnetic field. We explored the phase 
diagram of this system for a range of temperatures and 
magnetic fields, and studied the competition between the 
antiferromagnetic couplings and the external field; the 
crossover points where the magnetic field quenches the 
frustration are correctly captured, and the overall ex- 
perimental results are in good agreement with numerical 
simulations and theory. 
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Supplementary information: Digital Quantum Simulation of the Statistical Mechanics of a Frustrated 

Magnet 



Quantum circuit construction for creating the GETS where 



For clarity, we re-write the Hamiltonian of the frus- 
trated magnet as: 

H = J {siS2 + + S153) + (51 + S2 + S3) , (9) 

where s = {0, 1}, and Z \s) = ( — 1)^ \s). To construct the 
quantum circuit diagram for creating the GETS of the 
frustrated magnet, we first consider the general property 
for the following controlled operation: 

|5i52)0|O) ^ |siS2)0(cos6>|O)+sin(9|l)) , (10) 

where 

cos = yJe-^J^^^+^^)-(^^/K (11) 

and 



sin l9 = We+/5^(^i+^2)+/3/i/i^ 



and 



^ — g-^^(si+s2)- 



(12) 



(13) 



Note that the numerators in Eq. (11) and Eq. (12) are 
chosen such that they give the correct weight of the Boltz- 
mann factors. Now, one can show that 



(14) 



where 
A = 
b = 



2 Vcosh {2(3J^(3h) cosh {2f3J-f3h) cosh {(3h)^ 
cosh (2/3 J + Ph) cosh (2/3 J - ph) 



1/2 



4/3 



In 



cosh^ (Ph) 



1 /cosh(2/3J + /3/i) 
4p ^ \cosh{2pJ - ph) 



(15) 



Gombining these results, we can interpret the operation 
in Eq. (10) as the one that gives the correct Boltzmann 
factors to the terms involving the third qubit, but it will 



renormalize the first two qubits, according to Eq. (14). 

This suggests that in order to prepare the three-qubit 
GETS, we will need to prepare the GETS for the first two 
qubits with respect to a Hamiltonian which takes into 
account the renormalization effect caused by the third 
qubit: 



H = JsiS2 + (51 + S2) 



(16) 



and 



J = J-b 



h = h — c . 



(17) 
(18) 



This can be achieved by the gates Ux^ Uy and Uz shown 
in the quantum circuit diagram Eq. ([T])c. Their explicit 
forms are as follows: Ux = R{Ox)^ Uy = R{Oy), and 
Uz = R{ez)R{Oyy\ where 



and 



R{0) 



COS Or 



COS U — sm f 
sin 6 cos 6 



and 



cos 0,j 



and 





2 cosh 


'p {h-g)' 


g-/3(j+/i) 


2 cosh 





(19) 



(20) 



(21) 



COS Oz 



2 cosh 



Here 



In 



p(-J^h) 
cosh {pj + ph)j 



(22) 



(23) 



cosh {pJ - ph^ 
The second part of the quantum circuit diagram 



implements the transformation Eq. (10), and involves 
the gates Uq^ Ui and U2. Their explicit forms 
are given by: Uo = R{0o), Ui = ((9i) (l9o)"\ 
and U2 = R{02)T-^ {OuOo), where T {61,60) = 
R{6r)R{6oy'R{6i), 



cos ^0 



-i(32J+(3h) 



2 cosh {p2J ^ph) 



and 



and 



cos 61 



-(3h 



COS 62 



2 cosh (ph) 



I ^(32J-(3h 

2 cosh {p2J- ph) 



(24) 



(25) 



(26) 
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Partial recovery from isotropic noise 

We applied an empirical transformation of the mea- 
surement results to partially remove decoherence errors. 
There is no rigorous proof that this transformation will 
work for all circumstances, however, there are reasonable 
motivations for it. In this experiment, all of the observ- 
ables {Um) are obtained by measuring the coherence of 
the probe qubit (see Eq. Q). For the moment, if we only 
take into account the T2 decay of the probe qubit, then 
all of the measured values should be smaller than the ac- 
tual value by a value roughly equal to 7^ = e~^l^'^ . Here 
r is the evolution time, and we assume it is the same for 
all of the measurements. This assumption is equivalent 
to saying that the prepared GETS \^ ^) is subject to the 
depolarizing channel: 



^ ^ (Pc) = (1 - ^) ^ + W 



(27) 



where pc = 1^/3) ^ and I) = 8 is the dimension of the 
GETS. In the ideal case where 77 is known, the GETS can 
be perfectly recovered from by the following transfor- 
mation: 



/ 



V 



I 
D 



(28) 



Or equivalently, we may simply multiply all of the ob- 
servables by the factor r]. 

In reality, a depolarizing channel as a noise model is a 
pure assumption, and it is almost impossible to determine 
the exact value of 77, e.g. in the case where the decay 
rate for each observable may be different from the others. 
Nonetheless, partial recovery from the noise is possible 
when the decay rates are roughly the same, and we may 
estimate r] in an average sense. To see this, consider the 
trace-preserving transformation: 



D 



1 
A 



D 



V 

V 



V a; d 



(29) 



where A is assumed to be close to r]. It is easy to check 
that this transformation is not positive. However, as long 
as A is close to 77, the rescaled state p' is approximately 
equal to i-e., the original GETS. 

In fact, for the depolarizing channel applied to a pure 
state, i.e., Tr{pl) = 1, a systematic estimation for 77 is 
possible. To see this, consider taking the trace of the 
square of p^ in Eq. (27), we find that 



l Tr{pl)-l/D 
1 - l/D 



(30) 



When D is much larger than 1, we have the approxima- 
tion 



If we take this correction to the density matrix to cal- 
culate the fidelity with the \^ ^) GETS, the result, in 
the limit I) ^ 1, is the same as that of the projection 
between them: 



1^/3) 



(32) 



Although, in theory, Eq. (30) is an exact expression 



T] ^ ^JTr {pD . 



(31) 



for 7/, in practice, it can easily cause the final density 
matrix to have negative eigenvalues. Empirically, the 
approximation 77 ^ \/Tr (p^) works much better. In this 
experiment, we chose rj = 0.6316, which is equal to the 
square root of the purity of the measured GETS at /3 = 11 
and h = 5. With a single fitting parameter, the total 
magnetizations for all of the states are rescaled to values 
much closer to the theoretical values, as shown in Fig.[3^- 
d. 

The term "isotropic error" in the main text refers to 
the part of the error generated by the uniform part of 
the decoherence. To make the statement more quanti- 
tative, as an example, consider two observables which 
suffer from two different dephasing rates 71 and 72. The 
isotropic error can be quantified by defining the mean 
value jrn = (71 +72) /2. By rescaling the factor e^"^^ 
to both observables, the effective decay rate becomes the 
anisotropic error rates 71 — and 72 — 7m , which is zero 
when 7i =72. 

Lastly, we emphasize again that this decoherence 
model is a pure assumption, and to motivate the use of 
it, we considered only the dephasing of the probe qubit. 
In reality, there are also decoherence channels for the 
other qubits. On the other hand, the depolarizing noise 
model is a useful approximation when we consider that 
the applied pulse sequences may randomize the noise to 
some extent. In the limit that we apply a large number 
of random operations uniformly drawn from the Glifford 
gates, we would in fact get exactly a uniform depolariz- 
ing channel (See, e.g., G A Ryan et al 2009 New J. Phys. 
11 013034). 



Ti relaxation of the proton spins 

In implement ating the pseudopure state, ps starts from 
the state a^lllllim. 

In the experiment, the proton spins M, Hi and H2 
are prepared in the state OOa^. We measured the decay 
of the state OOa^MM in order to estimate the effect of 
the protons' Ti relaxation process on the measurement 
of the GETS. The state OOcr^llllllll is prepared from ps = 
a^lUm through phase cycling. After a delay time, a 
7r/2 readout pulse for H2 is applied. The intensity of the 
signal against the delay time is shown in Fig. |5| with a 
fit yielding a relaxation time of 2.95 s. 
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FIG. 5: The decay of the H2 signal obtained by a 7r/2 readout pulse to the state OOcr^llJlll. The experimental data are denoted 
by *, and the fitting result is shown as the solid curve. 



Supplementary data for various measurement observables, the spectra of the 

probe qubit, and selected measurement results not in- 
The supplementary figures are shown below. They in- eluded in the main text, 
elude a table summarizing the experimental durations 
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Measurement 


Experiment duration (s) 


Theoretical value 


Experimentally measured value 


{Zi) 


0.35 


-1/3 


-0.2497 




0.46 


-1/3 


-0.1208 


(^3> 


0.57 


-1/3 


-0.2564 


{ZxZ2) 


0.62 


-1/3 


-0.1335 


{Z2Z3) 


0.58 


-1/3 


-0.3377 


{Z1Z3) 


0.76 


-1/3 


-0.1549 




0.59 


1 


0.5110 


(Xi) 


0.37 





0.2007 


{X2} 


0.49 





0.0743 


(Xs) 


0.63 





0.0154 


(Yi) 


0.37 





-0.1261 


{Y2) 


0.49 





-0.1375 


(Ys) 


0.63 





-0.1769 



TABLE I: Experimental parameters for estimating errors in implementation. Experiment duration for measuring Um includes 
the preparation of GETS, which takes 0.33 s. The theoretical and experimental data correspond to T = 1/11 and /i = 1. The 
differences in the errors in for Um give rise to the anisotropic part of the error in the reconstructed density matrix, which 
cannot be removed by the empirical transformation described above. 




FIG. 6: Density matrices for T — 1/11 and h — 5. From left to right, the four columns show the results in theory, by simulation 
(including T2), in experiment by full state tomography [26^, and in experiment using the phase kick-back readout. The top and 
bottom rows show the real and imaginary parts, respectively. Gompared with the theoretical result, the state fidelities from 
the simulation, state tomography and phase kick-back method are 79.8%, 60.2% and 59.8%, respectively. 
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FIG. 7: (a) Spectra of Ci obtained by 7v/2 readout pulses when the system is in the labelled pseudo-pure state ps = OOcr^OOOO 
(red) and reference state cr^llcr^llllll (blue), respectively. The vertical axes have the same scaling (arbitrary units). (b,c) NMR 
spectra of Ci for measuring {Z1Z2Z3) in experiment (b) and by simulation that includes T2 effects (c), with /3 = 5, /i = — 1 
(red), (blue), and 1 (black). In each spectrum, only the peaks marked by M in state |0) are shown, because the peaks marked 
by M in state |1) are too weak to attribute observable signals), (d) For reference, simulated spectra of Ci marked by spin H2 
in state |0) (dashed) and |1) (solid) are shown (with C2-C4 in the state 1®^). 
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FIG. 8: Measurement results of (X), (Y) and (Z) for the high-temperature hmit /3 = 1 (left) and the low-temperature limit 
/3 = 11 (right). 



